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LONG-TERM  GOALS 

The  development  of  new  geoacoustic  inversion  methods,  their  use  in  the  analysis  of  shallow  water 
experimental  data,  and  evaluation  of  geoacoustic  model  and  parameter  uncertainties  including  the 
mapping  of  these  uncertainties  through  to  system  performance  uncertainties. 

OBJECTIVES 

The  development  of  new  geoacoustic  inversion  procedures  for  use  into  the  kHz  frequency  regime,  the 
development  of  methods  for  estimating  the  entire  posteriori  probability  densities  of  the  geoacoustic 
parameters  being  investigated  along  with  the  mapping  of  these  parameter  uncertainties  through  to 
characterizations  of  applied  interest  (e.g.  transmission  loss),  and  the  demonstration  of  their  use  in  the 
analysis  of  data  collected  during  the  Shallow  Water  2006  (SW06)  experiment. 

APPROACH 

Geoacoustic  inversion  involves  a  number  of  components:  (a)  representation  of  the  ocean  environment, 
(b)  the  inversion  procedure  selected  (e.g.  genetic  algorithm  or  simulated  annealing)  including  the 
forward  propagation  model  implemented,  and  (c)  the  estimation  of  uncertainties  associated  with  the 
parameter  estimates.  The  latter  is  critical  to  facilitate  the  mapping  of  these  uncertainties  into 
characterizations  of  applied  interest  including  the  prediction  of  total  system  performance. 

Substantial  experience  exists  in  the  application  of  full-field  geoacoustic  inversion  methods.  These  have 
been  implemented  in  a  number  of  geometries  (e.g.  fixed  vertical  and  horizontal  arrays,  towed  arrays, 
and  sonobuoys)  and  have  been  shown  to  work  well  at  low  frequencies  (<  1  kHz).  The  application  of 
these  methods  at  higher  frequencies  (into  the  few  kHz  frequency  regime)  is  at  an  early  stage.  New 
methods  are  required  which  are  robust  to  modest  geoacoustic  heterogeneity  (seafloor  parameters  as 
well  as  bathymetry)  and  temporal  fluctuations  (sound  speed  structure,  surface  waves,  and  array 
dynamics). 

The  reporting  of  geoacoustic  parameter  estimates  without  their  associated  uncertainties  is  of  limited 
value.  Of  substantial  greater  utility  is  the  complete  a  posteriori  probability  density  (in  general,  the  joint 
density  between  all  parameters  being  estimated).  One  significant  benefit  of  obtaining  accurate  a 
posteriori  densities  of  the  geoacoustic  parameters  is  the  potential  to  map  these  through  to 
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characterizations  of  applied  interest  (e.g.  transmission  loss)  in  order  to  quantify  those  uncertainties  as 
well. 

The  Shallow  Water  2006  experiment  took  plaee  in  July-September  2006  on  the  outer  edge  of  the  New 
Jersey  continental  shelf  in  approximately  80  m  deep  water.  Both  narrowband  and  broadband 
transmissions  (source  tows  and  stations)  were  made  over  a  wide  range  of  frequencies  (50  Hz  -  5  kHz) 
ineluding  detailed  measurements  of  seafloor  strueture  and  water  eolumn  variability.  These  data  are 
available  for  geoaeoustie  inversion  purposes  and  the  investigation  of  how  nuisance  parameter 
uneertainty  (e.g.  water  eolumn  sound  speed  variability)  eouples  into  seafloor  parameter  uncertainty. 

WORK  COMPLETED 

A  method  for  estimating  transmission  loss  (TL)  has  been  developed  whieh  ineorporates  uncertainties 
in  the  acoustie  environment.  Speeifieally,  an  approaeh  has  been  derived  and  validated  for  the  statistieal 
estimation  of  TL  based  on  the  posterior  probability  density  of  environmental  parameters  obtained  from 
the  geoaeoustie  inversion  process  [2,4]. 

Participation  in  the  Shallow  Water  2006  experiment  took  plaee  23  August  -  7  September  2006. 

Several  researeh  groups  were  on  board  the  RW  Knorr  (W.  Hodgkiss  and  P.  Gerstoft,  SIO;  D.  Knobles 
and  P.  Wilson,  ARL/UT;  R.  Chapman,  U.  Vietoria;  J.  Miller  and  G.  Potty,  URI).  Four  array  systems 
were  deployed  including  two  MPL  16-element  VLAs  and  two  ARL/UT  SWAMI  HLA/VLAs.  These 
array  systems  reeorded  a  variety  of  souree  tow  and  souree  station  transmissions  including  CW  tonal 
sets  (50-950  Hz  and  1.5-4. 5  kHz),  LFM  chirps  (100-900  Hz  and  1.5-4. 5  kHz),  light  bulbs,  and 
combustive  source  transmissions. 

RESULTS 

Previously,  geoaeoustie  inversion  results  from  the  ASIAEX  East  China  Sea  experiment  were  reported 
for  low  frequeney  source  tow  transmissions  (195,  295,  and  395  Hz)  [1].  Subsequently,  these  results 
were  used  in  an  initial  demonstration  of  an  approaeh  for  mapping  uncertainty  in  the  geoaeoustie 
parameter  estimates  into  uneertainty  in  predicted  transmission  loss  [2]. 

Quantifying  uneertainty  for  geoaeoustie  parameter  estimates  requires  estimation  of  the  uncertainties  in 
the  data  due  to  both  ambient  noise  as  well  as  modeling  errors  with  the  latter  aeeounting  for  simplistie 
assumptions  about  seafloor  strueture,  water  column  variability,  range  dependeneies,  ete.  [3].  Both  of 
these  eombine  in  a  total  error  varianee  whieh  describes  the  data  uncertainty. 

The  observed  data  corresponds  to  a  source  position  approximately  50  m  deep  and  1.7  km  away  from  a 
16  element  vertical  line  array  deployed  in  105  m  deep  water.  A  total  of  13  parameters  (geoaeoustie, 
geometric,  and  water  eolumn  sound  speed)  were  estimated  in  the  original  inversion  which  made  use  of 
a  genetie  algorithm  based  global  optimization  procedure  [1].  In  the  initial  approaeh,  exhaustive  grid 
sampling  was  used  to  obtain  the  geoaeoustie  uncertainties  and  map  these  to  the  TE  domain  [2].  This 
was  feasible  beeause  only  4  model  parameters  (instead  of  all  13)  were  explored.  However,  a  more 
realistie  inversion  will  have  a  larger  number  of  parameters.  Here  a  Markov  chain  Monte  Carlo 
(MCMC)  procedure  is  employed  in  the  inversion  process  to  sample  the  posterior  probability  density  of 
all  13  geoaeoustie  parameters.  Then  these  sampled  parameters  are  mapped  to  the  transmission  loss 
domain  where  a  full  multidimensional  probability  distribution  of  TL  as  a  funetion  of  range  and  depth  is 
obtained.  Based  on  the  geoaeoustie  inversion  results,  the  predicted  TL  and  its  variability  are  estimated 
and  then  eompared  with  the  measured  TL. 


Fig.  1  summarizes  the  estimation  of  TL  (usage  domain)  from  oeean  aeoustie  data  observed  on  a 
vertieal  or  horizontal  array  (data  domain).  The  geoaeoustie  inverse  problem  is  solved  as  an 
intermediate  step  to  obtain  the  posterior  distribution  of  environmental  parameters  p(m|d) 
(environmental  domain).  We  are  not  direetly  interested  in  the  environment  itself  but  rather  a  statistieal 
estimation  of  the  TL  field  (usage  domain).  Based  on  the  posterior  distribution  p(m|d),  the  probability 
distribution  of  transmission  loss  p(u|d)  is  obtained  via  Monte  Carlo  integration.  From  this  TL 
probability  distribution,  all  relevant  statisties  of  TL  ean  be  obtained  (e.g.  median,  pereentiles,  and 
eorrelation  eoeffieients). 

Using  this  proeedure.  Fig.  2  shows  the  posterior  distribution  of  TL  vs.  range  at  295  Hz  for  array  El  #7 
(69.5  m  depth).  Fig.  2(a)  shows  the  eontour  of  the  predietive  distribution  of  TL  vs.  range.  Gray  levels 
represent  the  probability  density.  Darker  shades  mean  higher  probability  of  observing  the  predieted  TL 
value.  Predietive  distributions  at  two  different  ranges  are  shown  in  Figs.  2(b)  and  2(o)  whieh 
eorrespond  to  regions  of  eonstruetive  and  destruetive  interferenee,  respeetively.  Sinee  the  distribution 
of  TL  is  often  poorly  approximated  by  a  normal  distribution,  the  eentral  tendeney  (median)  and  spread 
of  the  TL  distribution  (5**^  and  95*  pereentiles  or  90%  eredibibty  interval)  are  indieated.  Fig.  2(d) 
summarizes  the  predietive  distributions  by  the  median  (heavy  line)  and  the  90%  Cl  (gray  area).  This  is 
a  praetieal  way  to  eonvey  the  uneertainty  in  TL. 

These  predietive  distributions  are  eompared  with  aetual  TL  observations  in  Fig.  3  for  a  souree  depth  of 
approximately  50  m.  The  observed  TL  (dots)  are  eompared  with  the  predieted  TL  statisties  (solid  line 
with  gray  area)  for  the  frequeneies  195,  295,  and  395  Hz  (left  to  right)  and  for  array  Els  #1,  7,  and  16 
(bottom  to  top  eorresponding  to  depths  99.5  m,  69.5  m,  and  24.5  m,  respeetively).  In  general,  there  is 
good  agreement  with  the  pereentage  of  observed  TL  data  points  within  the  eredibibty  interval.  Thus, 
the  geoaeoustie  inversion  statisties  have  eaptured  most  of  the  uneertainty  in  the  environment. 

IMPACT  /  APPLICATIONS 

Geoaeoustie  inversion  teehniques  are  of  general  interest  for  the  estimation  of  waveguide  parameters 
thus  faeilitating  system  performanee  predietion  in  shallow  water.  Natural  transition  paths  for  these 
results  will  be  SPAWAR  (PMW-155)  and  NAVSEA  (ASTO). 

RELATED  PROJECTS 

This  projeet  is  one  of  several  sponsored  by  ONR  Code  3210A  to  partieipate  in  the  Shallow  Water 
2006  experiment  and  partieipate  in  the  analysis  of  the  resulting  data. 
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Data  domain  V 


Figure  1.  An  observation  d  is  mapped  into  a  distribution  of  environmental  parameters  m  that 
potentially  could  have  generated  it.  These  environmental  parameters  then  are  mapped  into  the 

usage  domain. 
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Figure  2.  Posterior  distribution  of  TL  vs.  range  at  295  Hz  for  array  El  #7  (69.5  m  depth):  (a) 
Contour  of  posterior  distribution  for  TL  vs.  range,  (b)  and  (c)  Posterior  distributions  of  TL  at  two 
different  ranges  (2. 75  km  and  2.85  km).  These  correspond  to  cuts  (vertical  dashed  lines)  through  the 
contour,  (d)  Statistics  of  the  predicted  TL  vs.  range.  The  solid  line  with  gray  area  around  shows  the 
median  and  the  90%  interval  of  the  posterior  distribution. 
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Figure  3.  Predicted  and  measured  TL  (dots)  for  array  Els  #1,  7,  and  16  for  frequencies  195,  295, 
and  395  Hz.  The  median  of  the  predicted  TL  (solid  line)  is  shown  together  with  the  90%  credibility 

interval  (gray  area). 


